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The longwall method of mining in underground coal seams is very efficient in uniform seams, but coal seam 
anomalies can make the method unprofitable and unsafe. This paper describes the theoretical basis for detection 
of coal seam anomalies using medium frequency (MF) radio transmission over paths on the order of 200 m in 
length. The key to the method is the sensitivity of the attenuation rate of the coal seam mode of propagation to 
changes in the coal seam parameters, such as height or electrical conductivity. From a large number of 
transmission paths, the principles of tomography can be used to reconstruct an image of the seam. 
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1. Introduction 

Most coal seams are horizontally bedded deposits on 
the order of a few meters in thickness, and an ideal seam 
is fairly flat with little variation in thickness over a large 
region. In the longwall method of coal mining [l] 1 , two 
horizontal parallel entries separated by approximately 
150 m are driven in the coal seam. The coal between the 
two entries is then mined out as the "longwall" retreats. 
The longwall method is very efficient for uniform coal 
seams, but coal seam anomalies can make the method 
unprofitable and unsafe. A remote sensing method that 
could detect coal seam anomalies would be extremely 
useful. 

Seismic methods are currently being explored in coal 
seams [2], but their effectiveness is as yet un- 
demonstrated. Underground radars using VHF (30 
MHz-300 MHz) have been used for some short-range 
applications, but they do not have sufficient range [3] to 
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probe the entire coal panel between the longwall entries. 
During the 1970s, it was found that communication be- 
tween loop antennas using MF (300 kHz— 3 MHz) was 
possible in coal seams for horizontal ranges up to several 
hundred meters [4,5]. The dominant mode of propaga- 
tion at MF in coal seams is TM (transverse magnetic) 
with the magnetic field horizontally polarized. Because 
the electric field is primarily vertically polarized with 
only a small longitudinal component, the mode is nearly 
TEM (transverse electromagnetic) and is commonly re- 
ferred to as the quasi-TEM mode or coal seam mode. 
And because the propagation constant of the quasi- 
TEM mode depends on the coal seam parameters (thick- 
ness and electrical properties), it should be a useful 
mode of propagation for remote probing of the coal 
seams. 

This paper's purpose is to explore the theoretical fea- 
sibility of using MF radio transmission in remote sensing 
of coal seams, and its remaining sections are organized 
as follows. In section 2, the propagation constant and 
the field distribution are studied for a uniform coal seam. 
The dependence of the attenuation rate on the seam 
parameters is of particular importance. In section 3, the 
excitation of the TEM mode by a vertical loop 
(horizontal magnetic dipole) is analyzed. In section 4, 



385 



transmission through a longwall coal panel is examined 
and the question of how to process the transmission data 
is explored. This is an area where further work is re- 
quired. In section 5, conclusions are presented, and areas 
for further work are suggested. 

2. Attenuation Rate and Field Distribution 
for the Quasi-TEM Mode 

The geometry for a uniform coal seam is shown in 
figure 1 . The seam has thickness 2h , and the coal permit- 
tivity and conductivity are e c and cr c respectively. The 
permittivity and conductivity of the surrounding rock 
are e r and cr r respectively. Free space permeability ju is 
assumed everywhere. 



(V 2 +k T 2 )H y =0, \z\>h. 



(2) 



The lowest order mode is even in z, and we also 
require that H y is propagating in the x direction and 
satisfies eq (2). Thus we can write H y in the coal 
(|z | <h) in the following form: 



H y =H Q exp(— j k c S x) cos (k c C z), 



(3) 



where S 2 +C 2 =l. 



H is an arbitrary constant, and S is a normalized propa- 
gation constant which must be determined from the 
mode equation. The fields are independent ofj>, and the 
x and z components of the electric field, E x and E y , are 
given by the following for \z \ <h : 
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Figure 1 — Geometry for a uniform coal seam of thickness 2 h . 

The time dependence is exp (ja)t), and it is suppressed 
throughout. The wavenumbers for the coal and rock, k c 
and k t , are given by 



& c =a)Vu, € cc and £ r =cj)Vu. e r 



(1) 



where e cc =e c —j cr c /o) 

and e rc = e r —j cr r /&). 

The lowest order mode is transverse magnetic (TM), 
and for propagation in the x direction it has only an H y 
magnetic field component. The mode equation for this 
lowest order mode has been given elsewhere [6,7], but 
we will derive it here because we need the resultant field 
expressions in the coal and in the rock. H y must satisfy 
the following Helmholtz equation: 

(7 2 +k 2 )H y =0, \z\<h 



E= ^]_lH, andEz=T J_LH x (4) 



ya>€ cc dz 



JOi€ cc dX 



Substituting eq (3) into eq (4), we obtain the following 
expressions for E x and E z : 

E x = —j C 7} c Ho exp {-j k c S x) sin (k c C z), 
E z = —S 7] c Ho exp {—jk c Sx) cos {k c C z), (5) 



where t] c = Vja /e cc . 

For z>h, H y must satisfy the Helmholtz equation in 
eq (2). Also we require that H y must decay as z goes to 
positive infinity. An appropriate form for H y for z >h is 



H y =A exp (—jk c Sx) exp (— j u z), 



(6) 



where u =Vk r 2 -k c 2 S 2 

and Im (u)<0. 

A is an unknown constant, and Im denotes imaginary 
part. E x and E 2 are again determined from Maxwell's 
curl equation: 

E _ -1 a^v 
uA 



o) e r 



and 



exp (— j k c S x) exp (—j u z) 



(7) 



Ez = T7Z — TTf = — T^Z ex P (-J kc S x) exp (— j u z). 



jcoe rc dx 



O) € r 
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Similar expressions for H y , E x , and E 2 can be obtained 
for z <h from the fact that H y and E z are even in z and 
E x is odd in z. 

Atz=h,H y and E x must be continuous which leads to 
the following pair of equations: ^. p,-|_ 



H Q cos (k c C h)=A exp (— y w A), 



a) e r 



From the first equation A can be written 



A = 



H cos (k c C h) 
exp (— y u h) 



(9) 



By substituting eq (9) into the second equation in eq (8), 
we obtain the following mode equation to be solved for 
C: 



j k c C tanh (j k c C h)-\-j u e cc /€ rc =0. 



(10) 



Equation (10) is consistent with earlier results [6-8], and 
is in a form which is easy to solve numerically. The 
solution of eq (10) by Newton's method [9] is detailed in 
the Appendix. 

The reason that coal seams support a quasi-TEM 
mode quite effectively is that the coal conductivity <r c is 
normally much smaller than the rock conductivity cr r 
[4,7]. The propagation constant T of the quasi-TEM 
mode is given by 



r=jk c s=jk c V\-c 2 



(ii) 



where Re (T)>0, 



C is determined from Appendix eq (A-3) and Re indi- 
cates the real part. The attenuation rate of the coal seam 
is the easiest quantity to obtain from transmission mea- 
surements, and is given by 8.686 Re (r) in dB/m. Addi- 
tional information is available from the phase velocity of 
the quasi-TEM mode, but phase measurements have not 
been attempted in coal seams. 

In figure 2 we show the attenuation rate of the quasi- 
TEM mode as a function of frequency for the following 
parameters: 

cr c =\0~ 4 S/m, €0/60=6, or r =10 -1 S/m, € r /£ =15, and 
2/z=2m. 

The permittivities are normalized to the free space value 
e . The computer code, written to generate the attenu- 
ation rate, was checked by comparing the results with 
those of Delogne [7], and in all cases agreed to graphical 
accuracy. As expected, the attenuation rate increases 
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Figure 2— Attenuation rate as a function of frequency. Parameters: 
cr c = 10- 4 S/m, ec/e =6, ^=10"' S/m, e r /e =15, 2h=2m. 

with frequency just as it does for a plane wave in an 
infinite rock medium. This indicates that the chosen 
frequency should be as low as possible to increase the 
range. However, there are two reasons for not choosing 
too low a frequency. For remote sensing, the resolution 
decreases as the frequency is decreased. Also, antenna 
efficiency decreases as the frequency is decreased. As a 
result of the tradeoff between attenuation rate and an- 
tenna efficiency, the maximum communication range 
between loop antennas in a coal seam was obtained for 
frequencies on the order of 500 kHz [5]. 

In figures 3-5, we show the attenuation rate as a func- 
tion of various seam parameters for a frequency of 500 
kHz. Since the goal is to detect changes in seam param- 
eters from measurements of attenuation rate, it is de- 
sirable for the attenuation rate to be sensitive to changes 
in the seam parameters. The most important thing to 
detect in practical applications is a decrease in seam 
thickness 2h. Figure 3 shows that the attenuation rate 
increases rapidly as thickness decreases below 1 m. This 
qualitative behavior can be predicted from Appendix eq 
(A-6). In figure 4 the attenuation rate is seen to increase 
as coal conductivity cr c increases for both a plane wave 
and the quasi-TEM mode. The behavior of the attenu- 
ation rate as a function of rock conductivity ar t is more 
complicated, as shown in figure 5. For large or r , the wall 
loss decreases as cr r is increased as indicated by the 
second term in Appendix eq (A-6). However, as cr T be- 
comes small, more of the energy propagates in the rock, 
and the attenuation rate decreases as ar r decreases. 
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Figure 3-5— From above, clockwise: Attenuation rate as a function of 
seam thickness. Parameters: f= 500 kHz, cr c =10 -4 S/m, e.c/eo—6, 
o- r =10 _I S/m, e r /eo=15. Next, attenuation rate for both a plane 
wave and the quasi-TEM mode as a function of coal conductivity. 
Parameters: /=500 kHz, e c /e =6, o- r =10 _1 S/m, £ r /e =15, 
2/j =2m. Finally, attenuation rate as a function of rock conductivity. 
Parameters: /=500 kHz, cr c =10 -4 S/m, e.c/eo—6, e r /eo=15, 
2h=2m. 
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In figures 6-8 the electric and magnetic field distribu- 
tions are shown for three different values of rock con- 
ductivity <r r . The normalization field E is the value of 
E z at z =0 and is given by 



E =T) C SH . 



(12) 



All field components decay exponentially in the rock 
(z>h), and the decay rate is most rapid for the largest 
value of cr r . Inside the coal seam, H y and E z are the 
dominant field components, and they are nearly con- 
stant in z. Thus a vertical electric dipole or a horizontal 



magnetic dipole could be used to excite (or receive) the 
mode, and they would be insensitive to the vertical pos- 
ition. The horizontal electric field is zero at the center of 
the seam and is fairly small throughout the seam. Thus 
the mode is quasi-TEM. However, in the rock \E X \ is 
larger than | E z | because E z drops discontinuously to a 
small value in the rock. Thus the electric current flow is 
primarily horizontal in the rock. In addition to provid- 
ing information on excitation of the quasi-TEM mode, 
the field distributions in figures 6-8 also indicate how 
anomalies in the seam will be illuminated. 
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Figure 6— Electric and magnetic field distributions of the quasi-TEM 
mode. Parameters: /=500 kHz, <r c = 10 -4 S/m, ec/eo=6, cr r =10 -1 
S/m, e r /e =15, 2h=2m. 

3. Excitation of the Coal Seam Mode 
by a Vertical Loop Antenna 

From the field distributions in figures 6-8, it is clear 
that either a vertical electric dipole or a horizontal mag- 
netic dipole will be effective in exciting the quasi-TEM 
mode. Short electric dipoles (or monopoles) located in 
or near a conducting medium are generally inefficient 
because of near-field losses [10-12], The near-field losses 
are generally smaller for magnetic dipoles or small loops 
[11,13], because the dominant near field is the quasi- 
static magnetic field. Consequently loop antennas have 
been considered best for communication in coal seams 
[4,5]. In this section we analyze the excitation of the 
quasi-TEM mode by a horizontal magnetic dipole (a 
small vertical loop with its axis horizontal). 



Figure 7 — Electric and magnetic field distributions for a higher rock 
conductivity. 

We consider a small loop of area A carrying a current 
/ as shown in figure 9. The loop is centered on the z axis 
at z =z , and the loop axis is in the y direction. Thus the 
source can be considered a j; -directed magnetic dipole 
of moment IA . This source radiates both a continuous 
spectrum of plane waves and a discrete spectrum of 
waveguide modes. For sufficiently large horizontal dis- 
tances, the continuous spectrum can be ignored because 
it corresponds to highly attenuated waves traveling in 
the highly conducting rock walls. Also all higher order 
waveguide modes are highly attenuated because they 
are all well below cutoff. This is because | k c 1 2h is much 
less than unity at MF. Thus only the quasi-TEM mode 
which has no low frequency cutoff is significant at MF. 

The excitation of the quasi-TEM mode has been con- 
sidered by Delogne [7], and his result for the azimuthal 
component of the magnetic field H$ is 

k 2 1 A 
H$= ° , A cos (k c C z) cos (k c C z ) K { (Tp) cos <f>, 
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Figure 8 — Field distribution for a lower rock conductivity. 
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where cos <J>=— p=Vx 2 +y 2 , 



(13) 



A = 



k r 2 -k c 2 sin(k c C2h) 
+ k 2 -k c 2 C 2 (k c C2h) 



K x is a modified Bessel function [14], and C is deter- 
mined from the solution of the mode eq (10). Equation 
(13) is valid when both the source and the field point are 
located in the coal seam (\z \ <h and \z Q \ <h), and this 
is the usual case of interest. When | T \ p is large, we can 
replace Ki(Tp) by its asymptotic expansion [14]: 

*i(rp)~V^exp(-r P ). (14) 

21 p 



Using eq (14), we can rewrite eq (13) as 



H^-H^f- V-?L_ A cos (k c C z) cos 
* h 2Tp 

(h e C z ) exp [-j k c (S - l)p], (15) 



where 



I A k 2 
#i= 4 c exp (-j k c p) cos (j>. 

The reason for the normalization in eq (15) is that Hi is 
the far field of a magnetic dipole in an infinite coal 
medium for the case z =z . Note that the cos <J> factor 
leads to maximum in the plane of the loop (y =0) and 
nulls along the loop axis (x =0). The remaining factors 
yield the effect of the surrounding rock walls. The first 
factors involving p and h indicate that the coal seam 
mode has cylindrical spreading whereas H x has spherical 
spreading. A relates to the excitation of the mode. The 
two cosine factors are the "height-gain" factors for the 
source and the observer within the seam. The ex- 
ponential factor arises because the propagation constant 
of the mode is not the same as that of the infinite coal 
medium. 

For the case where the rock conductivity is much 
greater than the coal conductivity ( | k t \ > | k c | ), the ex- 
pression in eq (15) simplifies considerably. In this case, C 
and S are given approximately by Appendix eqs (A-5) 
and (A-6): 



A/ITandS^l+TTTT. 
k t h J 2k * h 



(16) 



Thus C is small, and S is close to unity. We assume that 
k c C h is small compared to unity, and thus the cosine 
factors are nearly unity: 



Figure 9— Vertical source loop (horizontal magnetic dipole) in a coal 
seam. 



cos (k c C z)=*cos (k c C z )— 1- 



(17) 
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This means that the height of the transmitter and re- 
ceiver within the coal seam is not important, and this is 
consistent with the numerical results in figures 6-8. The 
exponential factor can be simplified from eq (16): 



exp [-j k c {S - l)p]«exp ( 2^ )• (18) 

The excess attenuation over that of an infinite coal me- 
dium is given by eq (18), and is seen to be inversely 
proportional to that coal seam thickness 2h . This is an 
important result in remote sensing of the seam height. 
The remaining factors in eq (15) can also be simplified 
from eq (16), and the result for H$ is 



7T 



* h 2jk c h 



exp 



\2k t h)' 



(19) 



4. Interpretation of the Transmission Data 

In longwall mining there are normally two parallel 
entries approximately 150 m apart, and these entries are 
accessible for radio transmission measurements. A log- 
ical measurement scheme is to step the transmitting and 
receiving loop antennas along the entries as shown in 
figure 10. If there are N transmitting and iV receiving 
positions, then there are N 2 transmission measurements. 
The loops are assumed to have their axes in the y direc- 
tion, and the received voltage v is proportional to H^ 



v=j oi ]i Q A t H x , 
where H x =Ha cos d> 



(20) 



and A t is the area of the receiving loop. 

There are many possible methods for interpreting the 
transmission data. Since the coal seam supports only a 
single mode at MF, the attenuation rate of the quasi- 
TEM mode is the most logical quantity to work with. If 
we assume that the coal seam parameters vary slowly, 
we can approximately neglect reflection and refraction 
and assume straight ray propagation as indicated in fig- 
ure 10. In that case the magnitude of H x for the / th ray 
path can be approximately written: 



Iflrfh ^'^fr exp [-/ a (x,y) dp], (21) 
Vp, o 



where H t is a constant depending on the transmitting 
antenna and the coal seam parameters, a(x, y) is the 
attenuation rate in nepers/m and depends on x and y, p, 
is the length of the / th ray, and <f>, is the angle of the / th ray 
with the x axis. The form of eq (21) is justified from the 
theory for the uniform coal seam in the previous section. 
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Figure 10— Top view of a longwall mine. The transmitting and re- 
ceiving loops have their axes in the y direction and are moved along 
the parallel entries. 

From eqs (20) and (21), the magnitude of the received 
voltage v, for the / th ray path can be written: 



^ Vocosfr _j a(pcy)dpl (22) 

Vp, o 



where V is a constant which depends on both antennas 
and the seam parameters. Taking the logarithm of both 
sides of eq (22), we have 



r / \ j 1 Pq cos2 4>/ 



(23) 



In geophysical tomography [15-17], the usual approach 
is to divide the intervening region into some large num- 
ber of cells over which a is assumed to be constant. In 
that case the path integral in eq (23) can be reduced to 
a sum, and eq (23) can be written 



2D il X 1 =Y i , 



(24) 



where y,=ln 



V cos 2 (f>, 

IvilVp/ 



Xi equals the unknown attenuation rate a for the 1 th cell, 
Du equals the length of the z' th ray through the 1 th cell as 
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shown in figure 11, and |v,| is the magnitude of the 
measured voltage for the / th ray path. Since eq (24) 
applies to all rays, the total system of equations can be 
written in the following matrix forms: 



[A/] ra = [7/]. 



(25) 



[Yt] is a column vector with iV 2 elements corresponding 
to the iV 2 transmission measurements in figure 10. The 
number of elements M in the unknown column matrix 
[Xi] equals the number of cells, and normally M is cho- 
sen equal to or somewhat less than N 2 . 



cross-borehole probing [15-17], the past work on geo- 
physical tomography is directly applicable. The two 
differences are that the coal seam mode experiences cy- 
lindrical spreading rather than spherical spreading and 
that the unknown is the attenuation rate of the coal seam 
mode rather than the bulk attenuation of the medium. 
The cylindrical spreading results in the p,~ 1/2 factor in 
eq (24) rather than p, _1 for spherical spreading. 

Since the solution of eq (25) for [X t ] yields the attenu- 
ation rate a for the / th cell, there is a further question of 
what seam parameters actually produced that value of 
a. The attenuation rate a is given by 



a=-lm(k c S). 



(27) 




In general S depends on the seam parameters in a 
i Ray complicated manner via the mode eq (10). Even when 
the rock walls are highly conducting, a is fairly compli- 
cated: 

a— Im fe)+Re (^). (28) 



Figure 11 — Geometry for the / th ray. Xj is the attenuation rate in the/ h 
cell. 



In theory [X t ] could be determined from eq (25) by 
matrix inversion or pseudo-inversion: 



[X l ] = [D il ]- l [Y i ],M=N 2 



or 



(26) 



[X,] = 



[A/f[A/] 



[D il ] T [Y i ],M>N 2 , 



where T indicates transpose. 

In practice, the inversion of the linear equation in eq 
(25) has several difficulties. The values of [Y f ] are not 
known exactly because |v,| contains noise and mea- 
surement error. Both the path integral model in eq (23) 
and the discrete version in eq (24) are approximate. 
Also, the system of equations in eq (25) may be very 
large. For example, if the antenna spacing is 5 m over a 
length of 150 m, then the number of equations is 
30x30=900. In spite of these difficulties, there are nu- 
merous techniques which have been successfully ap- 
plied to geophysical tomography problems [15-17]. 
Some of the popular methods and their acronyms are: 
algebraic reconstruction technique (ART) [18], simulta- 
neous iterative reconstruction technique (SIRT) [19], 
and back projection technique (BPT) [20]. The descrip- 
tion and application of these methods to geophysical 
tomography is well covered [16,17]. Since our system of 
equations for the coal seam geometry as described by 
eqs (24) and (25) is nearly identical to the equations for 



Since a in eq (28) depends on three parameters (k c , k T , 
and h), we cannot determine the three parameters from 
a. However, in some cases additional geologic informa- 
tion might be available. For example if the rock and coal 
properties are known and constant, the seam height can 
be obtained from: 



2h^ 



Re (k c /k T ) 
a-flm(£ c )" 



(29) 



If eq (28) is not valid, 2h can be obtained from a curve 
as in figure 3. 

Even if it were not possible to determine the coal 
seam parameters, just the information that seam attenu- 
ation a changes in a region would be valuable to miners. 
Such a change would indicate that the seam is not uni- 
form, and core drilling could be done in that area to 
provide further information. In fact, in many cases the 
full tomographic processing of the transmission data 
might be more complicated than necessary. In some 
cases, anomalous regions have been readily apparent 
from the measured data because the few rays which pass 
through the region have shown very high attenuation 
[21]. The anomalies have later been verified by drilling 
or mining into the anomalous region. High conductivity 
anomalies in copper mines have also been found by a 
similar ray analysis [22]. 

The conditions under which the straight ray tomog- 
raphic approximation is valid have been discussed by 
Dines and Lytle [16], and similar conditions apply here 
if we replace the wavelength by guide wavelength and 
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the skin depth by the 1/e distance for the mode. Then 
the conditions become [16]: 



tions. The same difficulties occur for other geophysical 
tomography applications when strong scattering is 
present in either electromagnetic or seismic cases [25]. 



and 



^>RelS^)' forall/ '' 
-2Im(<S^ c )<Re(S^ c ). 



(30) 



The first condition is well satisfied at 500 kHz for 
pi > 100 m. The second condition is only marginally sat- 
isfied at 500 kHz, but is better satisfied when either the 
frequency or the coal dielectric constant e c is increased. 
Even when eq (30) is satisfied, it is also necessary that 
the seam parameters vary slowly in order that the 
straight ray approximation is valid. Otherwise reflection 
and refraction occur, and they can be included, but only 
with significant increase in complexity [23]. When scat- 
tering becomes very important, even further changes in 
the processing are required [24], and in some cases it is 
useful to search for the scattering pattern of expected 
anomalies [25]. More experience with coal seam anoma- 
lies is necessary to determine the best method for pro- 
cessing data. There is also the possibility that phase in- 
formation could be useful. 



5. Conclusions 

The feasibility of using the quasi-TEM mode at MF 
for remote sensing of coal seams has been studied. The 
mode has sufficient range for transmission between 
longwall entries which are normally separated by about 
150 m. The attenuation rate is sensitive to the coal seam 
parameters (coal conductivity, rock conductivity, and 
seam thickness), as shown in figures 3-5. Thus the mode 
should be useful in sensing variations in these parame- 
ters. The mode can be excited efficiently by a vertical 
electric dipole or a vertical loop (horizontal magnetic 
dipole). 

Transmission between vertical loops located in the 
two entries appears to be a useful method for probing 
the coal seam. Because a coal seam supports only a 
single mode which is nearly TEM, the system of linear 
equations to be solved is nearly identical to the system of 
equations in geophysical tomography. Thus the work 
which has been done on geophysical tomography for 
cross-borehole probing is directly applicable to the 
probing of coal seams. In some cases, such as a fault or 
washout of the seam, the seam properties change rap- 
idly, and tomography is not valid; strong scattering and 
reflection occur, and some other type of processing is 
required. However, the transmission measurements will 
still show that the seam is not uniform, and this informa- 
tion might be sufficient for practical mining applica- 
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Appendix — Solution of Mode Equation 

To solve eq (10) by Newton's method, we first define the function /(C): 

/(C)=y k c C tanh (j k c C h)+fe) Vk 2 -k c 2 (l-C 2 ). (A-l) 

The derivative of /(C) with respect to C is 

/'(C)=^=y k c tanh (j k c C h)~k c 2 C h sech 2 (j k c C h) 



+j (€c/e r )*c 2 CV Vk 2 -k c 2 (1 -C 2 ). (A-2) 

If C„ is an n th estimate of the solution to /(C) =0, then by Newton's method [8] the n + 1 estimate is 

The iteration in eq (A-3) is continued until the change in C is sufficiently small. 

An initial estimate C can be obtained by considering the perfectly conducting parallel plate case 
(e rc = oo ). In this case, the second term in eq (A-l) is zero, and the solution is C =0. For large € rc /e cc , 
C is not zero, but it is still small. Thus we can replace the tanh function by its argument: 

j k c C (j k c C Q h) +j k„ CccA^O. (A-4) 

The solution to eq (A-3) is 

Co^VJl. (A-5) 

k T h 

Equation (A-5) is consistent with an earlier result from Wait [8], and we use it for our initial estimate 
in eq (A-3). Equation (A-5) can also be used to obtain an approximate propagation constant of the 
mode: 

1 



j k c S~j k c -V\ - C 2 ~y k c {\ -^ Co 2 ) 

■**. + 2ih- (A - 6) 

The first term j k c is the propagation constant of a plane wave in an infinite coal medium, and the 
second term is a correction to account for the finite conductivity of the surrounding rock. 
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